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Abstract. - The second hyperpolarizability 7jv(— 3cij; tj,a;,a;) of N double-bond finite chain 
of trans-polyacetylene is analyzed using the Su-Schrieffer-Heeger model to explain qualitative 
features of the size-dependence behavior of 7jv. Our study shows that ^n/N is nonmonotonic 
with N and that the nonmonotonicity is caused by the dominant contribution of the intra- 
band transition to 7jv in polyenes. Several important physical effects are discussed to reduce 
quantitative discrepancies between experimental and our results. 
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C The size saturation behavior of the second-order hyperpolarizabihty 7Ar[— (wi + 102 + 

^V ■ W3); oJi, UJ2, W3] of finite conjugated polymer (especially simple polyene chains) has been exten- 
sively studied both experimentally [1-5] and theoretically [6-21] for the past several decades. 
vQ ' A power law on molecular size with variable power exponent is often used to describe the 
f^ I magnitude of the off-resonant nonlinear response in scaling. That is, ^n ~ iV^, where N is 
the number of double bonds [2,3]. The scaling is expected to saturate for the large N (ther- 
2 ! modynamic limit), i.e., b = 1 as N ^ 00. Several measurements on cubic optical nonlinearity 
in long-chain polyene oligomers [2,3] have shown the following three main features of the size 
' , I dependence: (i) the exponent b ^ 3 — 3.5 for small N; (ii) 'Jn/N is nonmonotonic as a function 
(-H ■ oi N, that is, 'Jn /N increases with N initially, and after having reached a maximum value, 
Q I then gradually decreases (about 10%) to the saturation value as N increases; (iii) the onset 
O ■ of the saturation occurs at A^ ~ 90 by the nonmonotonic fitting. We remark here that the 
nonmonotonicity of 'Jn/N has been overlooked by most theoretical computations, and that 
in Rcf. [3] the experimental data were fitted by a monotonic curve within experimental error 
. ^ and the onset of the saturation occurs at iV ^ 60 by the monotonic fitting, 
j^ ' Theoretical studies range from simple tight-binding such as Hiickel or Su-Schrieffer-Heeger 

(SSH) models [6-11] to Parier-Parr-Pople (PPP) [14,15], Hubbard [16] and electron-hole pairs 
models [17,18]. Detailed quantitative comparisons between experiments and theories are quite 
difficult due to the following several reasons. First, most theoretical studies have been limited 
to planar irans-polyenes, while the polyene chains in the solution are usually 'disordered' 
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or 'worm-like' [3]. Second, experiments used series of polymers that, although related to 
polyenes, contain a central phenyle ring and incorporate rings that are saturated aside from 
a double bond in the "polyene" [3]. Third, the optical gaps are around 2.3 eV rather than 
typical polyacetylene's 1.8 eV in most theoretical studies [3]. Hence, even under the same 
assumption of the planar structure of polyene chains, different models coupled with various 
(analytic and/or numeric) approximations have predicted vastly different jn and saturation 
behaviors for polyenes [6-21]. 

In this letter, based on the SSH model and our previous work [22,23], we have derived an 
exact expression for the second-order hyperpolarizability of third-harmonic-generation (THG) 
7jv(— 3w; w, w, w) for the finite chain of N double bonds. By choosing the typical parameters 
for polyenes, we notice that the exponent b varies from 2 to 5 for small N (N < 21) in our 
computation. Our result also shows strong nonmonotonicity ofjj^/N versus N. These match 
the experimental observations (i) and (ii) mentioned above. Finally, we observe that there are 
large quantitative discrepancies between this theoretical study and experiments in observation 
(iii) and the magnitude of jn/N. 

Our computation shows that 7jv can be split into two parts: a positive part due to intra- 
band transitions and a negative part due to interband transitions. The positive contribution 
from intraband transitions dominates the negative contribution from interband transitions. 
A positive j^ in our computation is consistent with all existing polyene experiments [1-4] 
which show that -yjy is positive for any N. Moreover, we have carefully treated the so-called 
unphysical interference effects mentioned in the previous calculation [18]. Our results show 
that they can actually be identified as the boundary effects. Finally, since the zero- frequency 
limit 7Ar(0; 0, 0, 0)/A^ can not be measured directly, most experimentalists first measure the 
THG spectrum "/n{~3uj;uj,uj,uj)/N by choosing the frequency close to its three-photon res- 
onance [2,3] or the second-harmonic-generation spectrum 'yN{—'2uj;uj,io,0)/N [4], then apply 
the empirical extrapolation to obtain 7Ar(0; 0, 0,0)/A^. This empirical extrapolation largely 
depends on the three-photon-resonance of polyenes, but it also needs information on the linear 
absorption of solvent [2,3]. Here we have derived the exact expressions for both the static hy- 
perpolarizability (lo ~ 0) and dynamic hyperpolarizability THG {lu =^ 0). This could provide 
us some physical examinations of those empirical extrapolation methods. Despite the igno- 
rance of Goulomb interactions [1,12-18] and some other important effects (such as the effects 
of end groups [5] , the conformational disorder [24] , the solvent effect [25] , the segments or short 
conjugation length [20,21], etc) that must be considered in the more accurate quantitative 
calculations on one-dimensional(lD) polymers, our results show that the non-interacting SSH 
model (which is only first approximation of real physical systems) can nevertheless provide a 
clear understanding of the qualitative physical pictures of the saturation behavior. 

The Hamiltonian for the SSH model [26] is given by: 



H = -y2[to + (-1)'2H (AVi.sQ,. + CLCi+i.s) + 2nKu^ (1) 



Ls 



where to is the transfer integral or hopping between the nearest-neighbor sites, u is the dimer- 
ization displacement, a is the changing rate of the hopping, n is the total number of CH 
monomers, K is the elastic constant, and Cj ^(Ci^s) creates (annihilates) an tt electron at site 
I with spin s. Each site is occupied by one electron. With the lattice constant a and the 
definition of gap parameter A = Aau, we have the eigenenergies: 



e^{k) = -£c(fc) = -\/[2tocos{ka)f + [Asin{ka)f, (2) 

where £^(A:) and £c(fc) correspond to eigenenergies in the valence and conduction bands, re- 
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spectively. To avoid Jahn- Teller effects [14,27], we work with the nondegenerate ground 
states of Hiickel chain with n — Am + 2 sites. Therefore, for a polyene chain with N(= n/2) 
double-bonds, the pseudo-momentum vector k = 0, ±TT/Na, . . . , ±(2A^ — l)7r/(4A^a). 

We resolve this size dependence problem via the following two steps. First, for each N we 
fix ipj O! and K to find the minimum energy of the ground state in Eq. ^by varying u. When 
the minimum energy is achieved, we denote the corresponding value of uhy uq. Second, under 
the static lattice configuration determined by uq, we obtain an exact expression for j^ for 
the dimerized Hiickel chain with N double bonds. This two-step treatment takes into account 
the fact that n electrons make a much greater contribution to the 7jv in polyene chain than 
nuclei and a electrons [1]. The Peierls instability in the quasi-lD polymer system [28] is also 
carefully considered in this method. 

By choosing to = 2.5eV, a = A.leV/A, K = IS.Oey/A^ and a = 1.22A, we obtain the 
energy gap Eg ~ 1.89ey and mq ^ 0.058A when A^ — > oo. Both values are a little bit larger 
than those in the reported experiments where uq ^ 0.04A and Eg ^ X.'&eV [28]. However, the 
above parameters are not unreasonable since non-interacting models neglect the important 
Coulomb interactions in the ID polymer chains [28]. The inclusion of the strong Coulomb 
interactions are generally expected to reduce the dimerized constant Uq, hence it reduces the 
gap parameter A and the energy gap Eg [28] as well. 

TablelHshows the relationship between N , uq, 2A and Eg. We would like to point out that 
there is no dimerization distortion for the extremely short chain with A^ < 5. A Peierls- type 
distortion toward a bond-alternated geometric structure with 14 or more tt electrons has also 
been reported in the semiempirical calculations of cyanine dyes [1]. Since the dimerization 
distortion only exists for a reasonable size of polyene chains, we apply our calculation of 7jv 
only for N > 5. 

Following the same procedures that have been developed in our previous work [22,23] and 
by using the discrete summation of wave vectors k instead of integrals, we obtain the following 
expression for the THG hyperpolarizability for a single finite chain: 



" k{occ) 



1 



-(l-(5V)(a:2-l) 



216 252 24 12 



(3) 



X9(x2-(3z)2) a;7(x2 - (3z)2)2 a;9(a;2 _ ^2) 3.7(a;2_^2)2 



where S = A/2io, x = ec{k)/A and z = hu}/2A. Eq.^{2Na) recovers x^^H-3^;'^)'^)'^) for 
the infinite chain as A^ ^ oo [23] . The first two terms in Eq. O correspond to the contribu- 
tions from the interband transitions while the remaining terms correspond to those from the 
intraband transitions. 

For the finite chain, the expressions of x^"^^ or ^n/N also have another boundary term [29]. 
We have neglected this term in our study due to the following physical reasons. First, the 
boundary term disappears when A^ — > cxd. Second, it causes the so-called strong unphysical 
interference effects if it is included [18]. Third, the actual measurement of polyene materials 
in the solvent solution [2,3] are on chain but not on ring structures, and the random phase on 
two open ending groups should be expected. Therefore, the boundary term should not have 
any observable effect or play a strong role in real physical environments. 

The existing optical experiments are measured under the wavelength A = 1.91/xm (or 
Ulu ^ 0.65ey) [2,3]. This wavelength is very close to the edge of three-photon resonance. In 
order to guarantee the THG transition in the off-resonant region, we have chosen hco = 0.6eV 
in our calculation. Substituting the parameters in Table. HI into Eq. O we have computed 
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7Ar(— 3aj;aj,aj,aj)/A^ for hui = O.QeV and hu = OeV. The results are shown in Fig^jV and 
in Fig^3, respectively. On the same graphs, we have also compared them with the existing 
experimental results [3] . The graphs clearly show that the positive contribution from intraband 
transitions always dominates the negative contribution from interband transitions for any N. 
The fact that 7Ar is positive in our calculation agrees with the reported experiments [1]. It is 
also consistent with the theory presented by Agrawal et. al. [6] , but differs from Mclntyre et. 
al.^s [7] and Beratan et. aL's [8] results where only the negative interband contributions are 
considered. 

FigH shows jn/N firstly increases until N reaches about 20-30 and then decreases there- 
after. Though this trend is parallel to the nonmonotonic feature in existing experiments [2,3], 
it has never been reported in the existing theories [6-21]. ^n/N curve has a much longer 
tail and almost 2 orders larger value for hu) = Q.QeV than its static limit Tiio = O.QeV . It 
shows that the three-photon resonance plays a strong role to saturate at a relatively large 
N. 7Ar(— 3a;; a;,a;,a;)/A^ saturates at N2^^ ~ 60 for Tiu = O.GeV while at N]^f ~ 40 for its 
static counterpart. The much larger saturation length and jn value at finite freqencies than 
at the static limit has also been reported in Luo et. a/.'s work [30]. Our calculated N]^^ 
mismatches worse at OeV than at 0.6ey with the experiments that report N]^^ ~ 60 (without 
considering slight nonmonocity) for both cases [2,3]. The iV^j^ difference between this calcu- 
lation and experiment for OeV could be understood because the experiment uses the following 
extrapolation formula to obtain the static limit [2,3]: 



1 — i^max/^) 
1 ^ (SAmarr/A) 



7^(0;0,0,0)= ^ ;^^""'"^;/ 7^(-3^;^,^,^), (4) 



and A correspond to the wavelength of the maximum of the absorption in solution 
and the measurement, respectively. Eq.0] assumes that 7Ar(0; 0,0,0) is largely dependent on 
the three-photon process in 77v(— 3a;;a;, w, w). However, Eq. |3| shows that the one-photon 
process ((x^ — z^)~^ terms) in 7jv plays a much more important role at its static limit, even 
if it is very small compared to its three-photon counterpart at fno = 0.6el/. The one-photon's 
contribution was ignored in the above experimental extrapolation formula. 

The one-photon's influence on the saturation behavior in jn /N is similar to that of lin- 
ear susceptibility Xn ■ Fig[21is plotted according to Eq.(3.13) in Ref. [22]. Xn saturates at 
A^"^,^ ~ 25 ± 2. The non-monotonic jn/N up-and-down trend is dominated by the intraband 
contributions and are closely related with the linear saturation -/V^^. The intraband contri- 
butions are determined by gradient k terms or k-changing rate of dipole-transitions [6,23]. 
The intraband ^y™*'"'^ drastically increases with N until N reaches -/V^j, then the increasing 
of N (corresponding to the additional 7r-electron) lowers the average contribution of ^™*'''' 
until 7™*''''/A^ reaches its static limit. Meanwhile the negative but relatively weak interband 
rymter jncreases monotonically with N to its static limit. This may explain the up-and-down 
features oi'yN/N in experiments [2,3]. 

Our theoretical results of N]^^. ~ 40 and N"^^ ^ 25 at uj — are in agreement with Lu et. 
aL's ab initio calculations of polyene (iV^^j = 45 ± 5 and N^^^ = 25 ± 2) [19]. They are also 
comparable with the work of Shuai and Bredas {N]^f ~ 50) [9], Yu and Su (A^7ot ~ 50) [10] 
and Spano and Soos {N]^^. ~ 30 for 6 « 0.18) [11]. We have applied a factor of 1/2 on the 
number of n-site carbon to convert the results in Ref. [9-11]. 

FigOl shows the scaling power b versus N. b varies with both N and uj. It is also very 
sensitive to some other factors such as the hopping tg, gap parameter A, etc [9-11]. The 
scaling power law has also been extensively studied in many theoretical works [8-19]. Due to 
the fact that 6 is a derivative from jpf, here we only make a simple discussion. For N=9, we 
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have b = 4.83 at Q.QeV and b — 2.91 at the static hmit. The steadily decreasing trend in the 
graph is quite similar to Fig. 4 {5 — 0.18) in Ref. [11] when N < 30. 

Some quantitative discrepancies between this calculation and the existing experiments [2,3] 
are evident. Although a much large or upper bound limit of 7Ar value should be expected under 
the Hiickel model due to the ignorance of Coulomb interactions [11, 14], theoretical values of 
Jn/N are much larger than experimental data especially when huj — O.GeV. There are 
many factors to influence this result. One chief factor could be the conformational behavior 
in the solvents [3, 24, 25] and the fact that the polyene chain is no longer oriented in one 
direction. Hence the worm-like polymer chains may significantly reduce "fN value by averaging 
the contribution in 3-D space. The segment or short conjugation length treatment of polyene 
chains may also reduce 7jv [20,21]. Another important factor is the damping or the lifetime of 
excited states [9]. This may smear off the resonant peak and hence reduce jn- The optical gap 
A of experimental "polyene" is around 2.3 eV while ours is 1.9 eV. Since 'Jn/N ~ A~^ [6], 
theoretical values would be about (2.3/1.9)^ ~ 3.2 times as large as the experimental data due 
to this factor alone. Finally, the nondegenerate ground state resulting from the end groups or 
fixed ring structures [3-5] in experiments may also play a role. 

In conclusion, the study of the size dependence of 7jv based on the SSH model provides us 
a solid physical background to understand the saturation behaviors in polyene systems. Most 
valuable features of in the experiments of jn could be qualitatively explained under the schema 
of single-electron models. However, the quantitative comparison with the existing experiments 
in ID polymer system still shows that the SSH model is only first approximation for real 
physical systems. Further studies require the refinement of the model with the consideration 
of many other important factors such as the Coulomb interactions, conformational behaviors, 
segments effects of chains, solvent effects or the damping factor, etc. 
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Table I - Relationship between number of double bonds N , dimerized constant uo, gap parameter 
2A and actual energy gap Eg. (to = 2.5eV, a = A.leV/A and K = 18.0eV/A^). 

N uo(A) 2A(el/) Eg{eV) 
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0.0073 


0.238 


3.098 


7 


0.0481 


1.576 


2.704 
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0.0542 


1.779 


2.467 


11 


0.0562 


1.843 


2.314 


13 


0.0569 


1.868 


2.211 


15 


0.0572 


1.877 


2.140 


17 


0.0574 


1.882 


2.088 


19 


0.0574 


1.883 


2.051 


21 


0.0574 


1.884 


2.022 


23 


0.0575 


1.885 


2.000 


41 


0.0575 


1.885 


1.922 


81 


0.0575 


1.885 


1.894 


251 


0.0575 


1.885 


1.886 
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Fig. 1 - A. Computed THG 7iv(-3a;; cj,cj,cj)/iV at huj = 0.6e\/ (left Y scale in units of 10"^^ esu): 
total contribution (solid line), intraband contribution (long dashed line) and interband contribution 
(short dashed line) versus the experimental THG result jn{3uj)/N (Ref.3) at A = 1.91/im (solid square 
and the right Y scale in units of 10~''''esu); B. Computed static 7jv(0; 0, 0, 0)/A^ (left Y scale in units 
of 10"^'' esu): total contribution (solid line), intraband contribution (long dashed line) and interband 
contribution (short dashed line) versus the experimental extrapolation static 7jv(0)/A'^ (Ref.3) ( solid 
square and right Y scale in units of 10~^^esu). 
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Fig. 2 - Computed linear susceptibility x i^^^) versus A'^ at hcu = O.Oel^ and hu> — 0.6eV^ for 
N> 5. 
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Fig. 3 - Computed scaling power coefficient b{N,uj) = d[log'yM{—'iuj',i^,^,^)]/d[logN] versus logN at 
hu! = O.OeV (solid line) and hcj = O.GeV (long dashed line) for A'' > 9. 



